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1.  INTRODUCTION 

Shock-excited  lead  zirconute  titanate  (PZT)  ferroelectric  ceramic  materials  have 
been  used  in  single-shot  pulse  power  supplies  to  deliver  average  powers  that  range 
from  a few  hundred  kilowatts  to  a megawatt  in  a few  microseconds.1'3  The  bound 
surface  charge  is  released  from  the  poled  PZT  ceramic  as  a shock  wave  passes 
through  it.  Previous  experimental  investigations  of  charge  release  phenomena  under 
shock  wave  compression  have  been  concerned  primarily  with  short-circuit  conditions 
and  large  resistive  loads.2'4  A few  investigations  have  been  reported  for  capacitive 
loads1  and  inductive  loads.5  To  date  a detailed  theory  has  not  been  developed  for 
the  electrical  response  of  a shock-compressed  ferroelectric  ceramic  due  to  the 
complicated  nature  of  the  depoling  process.  The  purpose  of  this  report  is  to  derive 
electrical  expressions  for  the  response  of  shock-depoled  ferroelectric  ceramics  with 
various  loads  by  using  simple  assumptions  for  the  depoling  process.  Electrical 
expressions  are  obtained  for  the  charge,  voltage,  current,  and  energy  in  the  circuit 
for  resistive,  capacitive,  and  inductive  loads  that  range  from  short-circuit  to 
open-circuit  conditions.  Plots  are  presented  for  these  expressions.  The  results  will  be 
useful  in  estimating  the  actual  response  of  shock-depoled  PZT  materials  as  the 
magnitude  and  type  of  load  varies. 

PZT  ceramic  materials  are  usually  depoled  in  either  of  two  modes:  normal 
mode  or  axial  mode.  The  shock  wave  propagates  through  the  PZT  material  in  a 
direction  either  perpendicular  or  parallel  to  the  PZT  remanent  polarization  vector  lor 
normal-  or  axial-mode  depoling,  respectively.  Normal-mode  depoling  is  usually 
preferred  for  high-voltage  applications  since  for  axial-mode  depoling  time-dependent 
electric  fields  much  larger  than  the  dielectric  breakdown  strength  of  the  PZT 
material  can  be  generated  behind  the  shock  front.4  In  this  report,  only  normal-mode 
depoling  is  considered. 

The  simple  assumptions  that  are  made  concerning  the  shock  depoling  of  the 
PZ I material  and  the  two  PZT  geometries  that  are  considered  are  given  in  Section 
11.  Sections  III.  IV.  and  V contain  the  analyses  for  the  resistive,  capacitive,  and 
inductive  loads,  respectively.  Geometrical  size  estimates  of  PZT  units  for  pulsed 
power  supplies  for  selected  high-voltage  energy  converters  are  given  in  Section  VI. 
Section  VII  contains  the  summary. 
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II.  SIMPLIFYING  ASSUMPTIONS  FOR  THE  CHARGE  RELEASE 
FROM  A SHOCK-DEPOLED  PZT  UNIT 


* 


It  is  assumed  in  this  analysis  that  the  PZT  material  responds  in  what  will  be 
called  an  ideal  manner.  The  assumptions  for  this  type  of  response  are:  (1)  the 
shock  front  is  steady  and  moves  with  a single  shock  velocity  U,  (2)  the  permittivity 
of  the  PZT  material  is  strictly  constant  (no  dielectric  relaxation)  and  equal  ahead  of 
and  behind  the  shock  front,  (3)  no  electrical  breakdown  and  conduction  occur  in 
the  PZT  material,  (4)  the  PZT  material  is  depoled  completely  and  instantaneously  at 
the  shock  front  with  no  kinetic  effects  occurring,  (5)  the  strain  behind  the  shock 
front  is  negligible,  and  (6)  once  the  PZT  material  comes  under  shock-wave 
compression,  it  remains  in  that  stress  state.  Although  these  simplifying  assumptions 
do  not  hold  for  all  depoling  modes,  stress  values,  and  electrical  load  condiFons,4-6-7 
experimental  results  have  shown  that  a stress  range  exists  for  some  PZT  materials  in 
which  the  simplifying  assumptions  are  approximately  valid.1 ,8'10 

The  two  geometrical  shapes  for  the  PZT  material  that  are  considered  in  this 
study  are  the  rectangular  PZT  unit  and  the  hollow-cylinder  PZT  unit.  A schematic 
for  the  normal-mode  depoling  of  a rectangular  unit  is  shown  in  Figure  1(a).  A 
planar  shock  wave  is  produced  at  one  surface  of  the  unit  usually  by  explosive 
detonation  or  projectile  impact.  A schematic  for  the  normal-mode  depoling  of  a 
hollow-cylinder  unit  is  shown  in  Figure  1(b).  For  this  unit,  a shock  wave  is 
produced  by  axial  detonation  of  a concentric  cylinder  of  explosive  located  inside  the 
PZT  cylinder.  The  PZT  material  is  usually  surrounded  with  a high-dielectric-strength 
material  (such  as  alumina-filled  epoxy  or  transformer  oil)  to  reduce  the  possibility  of 
dielectric  breakdown  around  the  edges  of  the  unit  under  the  combined  effects  of 
high  voltage  and  stress.1-3 


For  the  rectangular  PZT  unit,  the  bound  surface  charge  is  released  in  a linear 
manner  by  the  propagating  shock  wave.  At  time  t the  released  charge  is 


\ P„A  j ■ OCKr. 

q(t)  =• 

/ P0A  , t > 7, 


where  A = wC  is  the  electroded  area  and  r = C/U  is  the  shock  transit  time.  The 
capacitances  behind  and  ahead  of  the  shock  front  are  given,  respectively,  by 


Ac  t 

C.  (t)  = . 0 < t < 7 . 

1 h 7 


1 


(b) 

Figure  1.  Schematic  ot  the  normal-mode  depoling  of  (a)  a rectangular  PZT  unit  and 
(b)  a hollow-cylinder  PZT  unit  with  an  arbitrary  load.  A unit  consists  of  a stack  of 
PZT  elements.  The  rectangular  unit  has  width  w,  length  t,  and  height  h.  A shock 
wave  enters  the  left  surface  of  the  unit  and  propagates  in  the  direction  C with  a 
steady  shock  velocity  U.  The  hollow-cylinder  unit  has  inner  radius  r0,  outer  radius 
r, , wall  thickness  and  height  h.  A shock  wave  enters  the  inner  surface  of  the 
unit  and  propagates  outward  in  a radial  direction  with  a steady  shock  velocity  U. 
At  time  t the  shock  front  is  indicated  by  the  dashed  lines.  Regions  1 and  2 are 
behind  and  ahead  of  the  shock  front,  respectively.  It  is  assumed  that  the  remanent 
polarization  P0  is  reduced  to  zero  by  the  propagating  shock  wave,  thereby  releasing 
the  bound  surface  charge  behind  the  shock  front.  The  released  charge  is  distributed 
between  the  PZT  unit  and  the  load,  resulting  in  a voltage  V(t). 

3 
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1 
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0 < t < r . 


(3) 


and 


(\<t) 


Ac 
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where  e is  the  constant  permittivity  of  the  PZT  material.  The  equivalent  capacitance 
of  the  PZT  unit  is  constant  and  given  by 


C = C.(t)  + C,(t)  = — , t > 0 . 

1 - li 

For  the  hollow-cylinder  PZT  unit,  the  shock  wave  propagates  outward 
radial  direction.  The  released  charge  at  time  t is  given  by 


(4) 

in  a 


q(t)  = 


0 < t < r . 


t > T . 


(5) 


where  A = nV.2 1 (2r(|/£)  + 1|  and  r = 2/U.  The  PZT  capacitances  behind  and  ahead 
of  the  shock  front  are,  respectively, 


and 


C,(t) 


Ac 

h 


0 < t < r . 


(6) 


(7) 
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The  equivalent  capacitance  of  the  PZT  unit  is  constant  and  given  by 


C = c,(t)  + c2(t)  = Y . t > 0 . (8) 

The  above  expressions  for  the  hollow-cylinder  PZT  unit  reduce  to  those  lor  the 
rectangular  PZT  unit  for  rQ/C  > 1.  In  this  case  the  geometrical  factor 
|(2r(t/8)  + t/T]/|(2r0/G)+ 1 1 in  Equations  (5).  (6),  and  (7)  approaches  a value  of  1. 
The  shock  depoling  of  a thin-walled  hollow-cylinder  unit  will  give  essentially  the 
same  time  dependence  of  charge  release  as  the  shock  depoling  of  a rectangular  unit. 


ill.  RESISTIVE  LOAD 

A.  GENERAL  CIRCUIT  EQUATIONS 


In  this  section,  the  load  in  Figures  1(a)  and  1(b)  is  taken  to  be  a resistance  R. 
The  voltage  balance  equation  is 


C 


i(t)R  = V(t) 


t > 0 , 


(9) 


where  qc(t)  is  the  charge  on  the  PZT  unit  at  time  t,  and  C is  the  capacitance  of 
the  unit  which  is  given  either  by  Equation  (4)  or  (8).  The  charge  conservation 
equation  is 

q(t)  = q . ( t ) + J i(t)dt  . t > 0 . (10) 

*o 

Differentiating  Equation  (10)  and  substituting  i(t)  from  Equation  (9)  gives 


dqc(t)  qc(t ) _ dq(  t ) 

dt  RC  dt 


(II) 
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This  equation  can  be  solved  for  q (t)  by  substituting  q(t)  from  either  Equation  (1) 
or  (5).  For  0 < t < t.  the  initial  condition  q .(0)  = 0 is  used  for  obtaining  qc(t). 
For  t > r,  the  derived  expression  for  qc(f)  is  used  as  an  initial  condition. 


B.  RECTANGULAR  PZT  UNIT 

In  this  section,  expressions  for  the  charge,  voltage,  current,  and  energy  for  the 
rectangular  PZT  unit  with  a resistive  load  are  derived.  The  charge  q.(t)  is  obtained 
by  solving  Equation  (11)  and  using  q(t)  from  Equation  (I).  The  result  is 


The  voltage  on  the  P7T  unit  can  be  obtained  from  the  expression 

q.  ( t ) 

V(t)  = -^r-  . t > 0 . (13) 

Figure  2(a)  is  a plot  of  the  normalized  charge  and  voltage  for  various  values  of 
RC/t.  The  circuit  current  is  obtained  from  the  expression 

q,.(t) 

i(t)  = , t > 0.  (14) 

Figure  2(b)  is  a plot  of  the  normalized  current  for  various  values  of  RC/t.  For  a 
given  shock  transit  time  t and  PZT  capacitance  C.  the  magnitude  of  the 
corresponding  resistive  load  R can  be  obtained  for  each  curve  in  the  figure.  Current 
histories  for  normal-mode-depoled  PZT  95/5  rectangular  units  have  been  measured  as 
a function  of  load  resistance.3  1 1 (PZT  95/5  is  a solid  solution  containing  95  mole 
% lead  zirconate  and  5 mole  % lead  titanate  with  niobium  as  a minor  added 
constituent.)  Although  charge  conduction  occurred  in  the  PZT  material,  the  current 
histories  are  similar  to  those  in  Figure  2(b). 
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NORMALIZED  CHARGE  OR 


TIME 

(a) 


TIME 

(b) 

Figure  2.  Electrical  relations  for  a rectangular  PZT  unit  with  a resistive  load  for 
various  values  of  RC/r.  (a)  Normalized  charge  qc(t)/PQA  or  voltage  V(t)/(PQ A/C)  on 
the  PZT  unit,  (b)  Normalized  circuit  current  i(t)/(PQ A/r). 
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The  amount  of  charge  that  has  (lowed  through  the  resistance  at  time  t is  given 
by  Equation  (10).  Making  substitutions  for  q(t)  and  qc(t)  gives  for  this  charge 


For  the  open-circuit  case  (RC/t  = °°),  q (t)  = q(t).  All  the  released  charge 
remains  on  the  PZT  unit  for  all  time  t.  For  this  case  the  charge  and  voltage  on  the 
PZT  unit  are  linear  functions  of  time  for  0 < t < t and  constant  for  t > r.  The 
current  is  zero  for  all  time  t.  Figure  2 shows  these  relationships.  For  the 

short-circuit  case  (RC/r  = 0),  qe(t)  = 0 for  all  t and  i(t)  = P()A/t  for  0 < t < r 
and  equals  zero  for  t > t.  No  charge  remains  on  the  PZT  unit.  All  the  charge  is 
released  in  the  time  interval  0 < t < r.  For  RC/t  between  these  two  extremes,  the 
PZT  unit  is  charged  during  the  interval  0 < t < r and  discharged  for  t > r.  As 

RC/t  increases  from  the  short-circuit  case,  more  charge  remains  on  the  PZT  unit 

during  the  shock  transit  time.  A maximum  voltage  of  P()A/C  occurs  at  shock  transit 
time  for  the  open-circuit  case.  For  a sufficiently  large  remanent  polarization  P(),  this 
maximum  voltage  can  be  large  enough  to  cause  electrical  breakdown  in  the  PZT 

unit.  The  area  under  the  current  pulse  curves  is  given  by  P(|A  for  all  values  of 
RC/t  except  °°.  This  means  that  for  large  enough  times,  all  the  released  charge 
(lows  through  the  resistance.  For  RC/t  = no  charge  (lows  through  the  resistance. 

It  is  not  necessary  to  have  exactly  zero  resistance  R to  observe  the  short-circuit 
current  pulse  in  Figure  2(b).  In  some  cases,  an  R of  100  £2  or  even  1 k£2  will  not 
affect  (as  observed  experimentally)  the  amplitude  of  the  short-circuit  current  pulse. 
It  is  only  necessary  that  RC/t  be  small  enough  that  the  shape  of  the  measured 
pulse  is  essentially  the  same  as  the  short-circuit  case.  For  example,  if  C = 250  pF 
and  t - 2.5  ps  (values  from  Reference  1)  the  ratio  RC/t  is  0.01  for  an  R of 

100  £2.  In  Figure  2(b),  the  current  pulse  for  this  case  would  be  essentially 

unchanged  from  the  RC/t  = 0 case. 

From  this  discussion  it  can  be  concluded  that  the  important  circuit  parameters 

are  the  shock  transit  time  t and  the  RC  time  constant,  t is  the  controlling 

parameter  for  small  RC/t;  charge  Hows  to  the  external  circuit  as  fast  as  it  is 
released  by  the  shock  wave.  RC  is  the  controlling  parameter  for  large  RC/t;  charge 
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is  released  to  the  PZT  unit  by  the  shock  wave  at  a faster  rate  than  it  (lows  to  the 
external  circuit.  Therefore,  charge  and  voltage  on  the  unit  increase  during  shock 
transit  time. 


An  interesting  point  to  consider  is  an  increase  in  the  magnitude  of  the  shock 
velocity.  As  the  shock  velocity  U approaches  infinity,  the  shock  transit  time  t 
approaches  zero.  This  means  that  there  is  no  charging  phase  for  the  PZT  unit:  i.e.. 
at  t = 0 the  charge  P()A  appears  on  the  PZT  unit.  Expressions  for  the  charge  q (t) 
and  current  i(t)  for  the  discharging  phase,  t > 0.  can  be  obtained  from 
Equations  (12)  and  (14),  respectively,  using  the  equations  for  t > r. 

For  t = 0.  the  charge  becomes 


q . ( t ,t  = 0 ) 


limo  P„A0^  W/RC  - l)e",/R(  = P()Ae',/R(  . t>0.  <l0) 


The  current  is  given  by 


i(t.T  = 0)  = 


JilV-t/RC 

RC  ^ 


t > 0 


These  equations  are  the  expressions  for  discharging  a capacitor  with  initial  charge 
P„A  through  a resistance  R.  For  R = q.(t.r  = 0)  = P(|A  and  i(t.r  = 0)  = 0; 
i.e..  the  charge  remains  on  the  PZT  unit.  For  R = 0,  q (t,r  = 0)  = 0 and 
i(t,T  = 0)  = «>.  This  case  is  analagous  to  short  circuiting  a charged  capacitor. 

An  expression  for  the  electrical  energy  dissipated  in  the  resistance  is  given  by 


(t)  dt  , t > 0 . 


To  solve  this  equation,  use  i(t)  = q (t)/RC  and  i(t)dt  = dq  - dq  . Substituting  gives 

i:,(t)=  J |rdq-  J |rdqc  . t>0.  (19) 


fMH 


Integrating  the  last  term  and  using  ET(t)  = E(.(t)+EL(t)  and 


Ec(t)  = ^ , t > 0, 


Ey(t)  = J ^ ’ X>  0 


This  is  the  electrical  energy  equation  for  the  circuit.  The  electrical  energy  in  the 
circuit  increases  for  0 < t < t to  a maximum  value  determined  by  the  integral  in 
the  above  equation.  For  t > r the  energy  does  not  increase  because  dq  = 0. 
Evaluating  the  integral  gives 


(q(t) - qc(t)|  , 


( Pn  A / RC\ 

) — jlq(t)-qc(t)l 

1 "r<t>  " ')  P„A  / RC  \ 

f — J[q(r)-  qc(r)l 


0 < t < T 


t > T . 


Figures  3(a),  3(b),  and  3(c)  are,  respectively,  the  normalized  energy  on  the  PZT  unit 
obtained  from  Equation  (20).  the  normalized  energy  dissipated  in  the  resistance 
obtained  from  Equation  (18),  and  the  total  normalized  circuit  energy  obtained  from 
Equation  (22). 

For  the  short-circuit  case  (RC/r  = 0).  Ej(t)  = 0 lor  all  t.  All  released  charge 
Hows  immediately  in  the  external  circuit.  E(. ( t ) = 0,  because  there  is  no  charge  or 
voltage  buildup  on  the  PZT  unit  and  EL(t)  = 0 because  the  resistance  is  zero.  This 
corresponds  to  the  case  of  current  flowing  (magnitude  P(|A/r),  but  no  electrical 
energy  is  being  stored  or  dissipated. 

Maximum  electrical  energy  is  stored  on  the  PZT  unit  for  the  open-circuit  case 
(RC/r  = °°).  Also,  EL(t)  - 0 for  all  t because  all  the  charge  is  stored  on  the 
capacitance.  As  RC/r  -*■  Equation  (22)  becomes 


(PnA)2  /t\2 


ET(t.RC/r  = oo) 


) (P„A)2 
' 2C 


0 ^ t < r 


t > r 
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Figure  3.  Energy  relations  for  a rectangular  PZT  unit  with  a resistive  load  for 
various  values  of  R Ch.  (a)  Normalized  energy  Ec  (t)/((PQ  A)2/2C)  on  the  PZT  unit, 
(b)  Normalized  energy  EL  (t)/((PQ  A)2/2C)  dissipated  in  the  resistive  load,  (c)  Total 
normalized  energy  ET  (t)/((PQ  A)2/2C). 


This  can  be  shown  by  using 


’r(lj)(r1)  +l'ish«r  order (^;)  . 


terms 


0 < t < r , (24) 


as  RC/t  -*■  °°. 


Tlie  expression  (P()A)2/2C  is  the  maximum  energy  for  an  ideal  PZT  unit.  This 
energy  is  usually  not  achieved  in  a depoling  experiment  for  a PZT  unit  poled  to 
maximum  remanent  polarization,  since  the  associated  voltage  P()A/C  gives  an  electric 
field  that  can  be  much  larger  than  the  breakdown  strength  of  the  PZT  material.  For 
example,  in  Reference  1,  some  of  the  PZT  56/44  units  that  were  shock  loaded  were 
cubes  with  a side  d of  12.7  mm.  The  units  had  a polarization  PQ  of  0.33  /aC/mm2 , 

an  available  charge  P A of  50  /aC,  and  a capacitance  C of  about  250  pF.  (PZT 

56/44  is  a solid  solution  containg  56  mole  % lead  zirconate  and  44  mole  7c  lead 

titanate  with  minor  added  constituents  of  niobium,  strontium,  and  lanthanum.)  For 

these  units  the  calculated  voltage  P(|A/C  is  200  kV  and  the  calculated  electric  field 


(P(|A/C)/d  is  15.7  kV/mm.  This  electric  field  is  about  four  times  larger  than  a 


breakdown  field  of  3.7  kV/mm  for  a similar  material,  PZT  65/35,  measured  under 
ambient  conditions.1-  (PZT  65/35  is  a solid  solution  containing  65  mole  % lead 
zirconate  and  35  mole  % lead  titanate  with  niobium  as  a minor  added  constituent.) 
To  prevent  electrical  breakdown,  it  becomes  necessary  to  choose  the  load  parameters 
so  that  the  circuit  voltage  and  the  corresponding  electric  field  are  reduced. 
Therefore,  only  a fraction  of  the  energy  (P0A)-/2C  and  energy  density 
((P0 A)2/2C)/d3  of  5 J and  2.5  J/cm3,  respectively,  are  available  from  the  PZT  unit. 
It  should  be  noted  that  if  the  PZT  56/44  material  were  poled  to  a remanent 
polarization  value  less  than  the  representative  value  of  0.33  gC/ninr . then  electrical 
breakdown  in  the  unit  would  be  less  likely  to  occur  and  the  energy  (P(|A)2/2C 
would  then  be  achieved  at  a lower  but  more  realistic  value. 


Figure  3 indicates  that  electrical  energy  is  both  supplied  to  the  PZT  unit  and 
dissipated  in  the  resistance  during  shock  transit  time  r.  But  for  t > r the  total 
circuit  energy  is  constant;  energy  on  the  PZT  unit  is  being  dissipated  in  the 
resistance.  For  example,  at  t = r for  the  RC/r  = 1 curves,  the  normalized  circuit 
energy  is  about  73 '/< . About  40%  is  stored  on  the  PZT  unit  and  about  33%  has 
been  dissipated  in  the  resistance.  For  t > r the  40%  stored  on  the  unit  is 
dissipated  in  about  2r.  Figure  4 is  a set  of  normalized  energy  curves  showing  this 
behavior  at  t = r for  RC/r  ranging  from  I0~3  to  103.  The  curve  for  the  energy 
dissipated  in  the  resistance  reaches  a maximum  for  RC/r  = 0.529.  The  maximum 
energy  dissipated  is  0.381((Po A)2/2C).  This  can  be  contrasted  with  the  maximum 
voltage  which  occurs  for  RC/r  = 00  and  the  maximum  charge  transfer  which  occurs 
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for  R C/t  = 0.  At  maximum  energy  transfer  the  energy  on  the  PZT  unit  is  about 
0.2((Po A)’/2C).  Also,  the  energy  stored  on  the  PZT  unit  equals  the  energy 
dissipated  in  the  resistance  for  RC/r  = 0.868. 


C.  HOLLOW-CYLINDER  PZT  UNIT 

In  this  section  expressions  are  derived  for  the  charge,  voltage,  current,  and 
energy  in  a circuit  consisting  of  a hollow-cylinder  PZT  unit  with  a resistive  load. 
Due  to  the  cylindrical  geometry,  the  expressions  are  more  complicated  than  those 
derived  in  the  last  section  for  the  rectangular  unit.  The  charge  q.(t)  on  the  PZT 
unit  is  obtained  by  solving  Equation  (11)  with  q(t)  from  Equation  (5).  This  gives 


This  equation  can  be  compared  with  Equation  (12|  for  the  rectangular  PZT  unit, 
l or  the  open-circuit  case  (RC'/r  = °°).  Equation  (25)  reduces  to 


This  is  just  Equation  (5)  for  the  charge  release  q(t). 


Equations  for  the  voltage  V(t)  on  the  hollow-cylinder  PZT  unit  and  current  i(t) 
in  the  circuit  can  be  obtained  by  substituting  Equation  (25)  for  q (t)  into 
Equations  (13)  and  (14).  respectively,  figure  5(a)  is  a plot  of  the  normalized  charge 
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Figure  4.  Normalized  energy  relations  for  a rectangular  PZT  unit  with  a resistive 
load  at  shock  transit  time  r as  a function  of  RC/r. 


and  voltage  on  the  PZT  unit  for  a representative  r()/C  value  of  1.  The  time 
dependences  for  the  cylindrical  and  rectangular  geometries  are  the  same  for  t > r. 
This  is  expected  since  after  the  shock  wave  passes  out  of  the  unit,  the  charge 
remaining  on  the  unit  in  both  cases  decays  with  the  RC  time  constant.  Notice  also 
that  the  curves  in  Figure  5(a)  are  similar  to  those  in  Figure  2(a)  for  the  rectangular 
unit.  In  addition,  it  can  be  shown  that  Fquation  (25)  reduces  to  Equation  (12)  for 
r() /V!  1.  This  can  be  shown  using  | ( 2 r0  / C ) - 2RC/r]  /[ ( 2r0 /C ) + 1 1 -*■  1 and 

l/|(2r0/C) - 2RC/r|  - 0 as  r0/C  -*  <*>. 

Figure  5(b)  is  a plot  of  the  normalized  current.  For  0 < t < r the 
short-circuit  current  increases  linearly  with  time.  This  increase  becomes  less  as  r ()/E 
becomes  larger  and.  as  expected,  approaches  the  case  for  the  rectangular  unit  as 


B T ) 


i(t,R  C/r  = 0)  = 


T "S  + 1 


0 < t < r 


t > t . 


At  t = r/2,  i(r/2,RC/r  = 0)  = P()A/r  as  shown  in  Figure  5(b).  The  integral  of  the 
short-circuit  current  is 


y i(t.RC/r  = 0)dt  = P„A 


as  expected. 

The  total  circuit  energy  can  be  derived  from  Equation  (21)  by  substituting  q(t) 
from  Equation  (5)  and  qc (t)  from  Equation  (25)  and  integrating.  The  result  is 


'"A(K<f T |q(t)-qc(t)| 

\t") 


2\ 

- q,.(0 


'r  V \A  3r  / />V 

— +4PnA  V - — I - ) . 0 

(!MC)  OX^H 
(;•:)  i 


P(i  A /RC’\\  vr  r /|  2q.(r)  \2i'  3/ 

( ) |q(r)  - (|  .<r)|  “ - * 4P«A  ~ — " 

cU(?-H  (?*-) 


It  can  be  seen  that  as  r()/8  -*■  00  this  equation  reduces  to  Equation  (22)  tor  the 
rectangular  unit.  For  the  open-circuit  case  (R C/t  = °°),  Equation  (29)  reduces  to 


(P, 


»a)2  («"  *i)  /t  v . 


0 < t < r 


2C 


Ex(t)  = 


-rn 


(30) 


(p„A)- 


2C 


t > r . 


To  show  this,  it  is  necessary  to  expand  qc(t)  in  Equation  (25)  for  RC/r 
then  form  the  expansion  q ( t ) - qc(t).  The  expansions  for  0 < t < r are 


00  and 


q,.U)  = 


PnA 


2ro 


, i) i . (m. „ ±\ Ji  . ... 

\e  t/t  \«  3t/  RCr  \3C  12r/(RC)-r 


0 < t < r, 

(31) 


and 


q(t) - qc(t) 


P..A 


2ro 


/r°  + ‘ \ t3  + ...1 

l\C  37/  RCr 

\3V  12r/  (RC)2r 

0 < t < r 


(32) 


Figures  6(a).  6(b),  and  6(c)  are.  respectively,  plots  of  the  normalized  PZT  energy, 
normalized  energy  dissipated  in  the  resistive  load,  and  total  normalized  energy. 
Figure  7 gives  the  normalized  energy  relations  for  the  hollow-cylinder  unit.  Plots  of 
total  energy,  capacitive  energy,  and  resistive  energy  are  given  at  t = r.  The  resistive 
energy  peaks  around  R C/t  = 0.5.  These  curves  are  similar  to  the  curves  in  Figure  4 
for  the  rectangular  PZT  unit. 
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It  can  be  seen  that  as  r()/C  -*  °°  this  equation  reduces  to  Equation  (22)  for  the 
rectangular  unit.  For  the  open-circuit  case  (RC/r  = °°),  Equation  (29)  reduces  to 


2r„  tV 


/ (P0A)2  \J tJ_  /t\2  , 0 < t < r 

2C  ° + iV  \ 

E T(‘)s  V K 7 

I tP  At2 


t > r 


To  show  this,  it  is  necessary  to  expand  qc(t)  in  Equation  (25)  for  RC/r 
then  form  the  expansion  q(t)  — q,(t).  The  expansions  for  0 < t < r are 


°°  and 


r/r  V 2 3r/  RCr  \3fi  12r/(RC)2r 


0 < t <r. 


q(t) ~ qc(t)  = 


^ + , 


\«  3r/  RCr  V 3V  12r/(RC)27 


0 < t < r 


Figures  6(a),  6(b),  and  6(c)  are,  respectively,  plots  of  the  normalized  IV.T  energy, 
normalized  energy  dissipated  in  the  resistive  load,  and  total  normalized  energy. 
Figure  7 gives  the  normalized  energy  relations  for  the  hollow-cylinder  unit.  Plots  of 
total  energy,  capacitive  energy,  and  resistive  energy  are  given  at  t = r.  The  resistive 
energy  peaks  around  RC/r  = 0.5.  These  curves  are  similar  to  the  curves  in  Figure  4 
for  the  rectangular  IYT  unit. 
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Figure  6.  Energy  relations  for  a hollow-cylinder  PZT  unit  with  a resistive  load  for 
various  values  of  BC/r.  Curves  are  given  for  rQ/V  = 1.  (a)  Normalized  energy 
Ec(t)/((P0A)2/2C)  on  the  PZT  unit,  (b)  Normalized  energy  EL  v *:)/(( PQ  A)2  /2C) 
dissipated  in  the  resistive  load,  (c)  Total  normalized  energy  ET  (t)/( (PQ  A'2 /2C) . 
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Figure  7.  Normalized  energy  relations  for  a hollow-cylinder  PZT  unit  with  a 
resistive  load  at  shock  transit  time  r as  a function  of  RC It.  Curves  are  given  for 

ro/4!  = 1- 


IV.  CAPACITIVE  LOAD 


A.  GENERAL  CIRCUIT  EQUATIONS 

In  this  section,  the  load  in  Figures  1(a)  and  1(h)  is  a 
voltage  balance  equation  is 


capacitance  CL-  The 


Melt)  = qL(t) 

C ' C, 


t > 0 


where  C is  the  capacitance  of  the  PZT  unit  given  by  either  Equation  (4)  or  (8). 
q{.  (t)  is  the  charge  on  the  PZT  unit  and  q(  (t)  is  the  charge  on  the  load  capacitor. 
The  charge  conservation  equation  is 


q(t)  = qc(t)  + qL(t)  , t > 0 , 


(34) 


where  q(t)  is  the  charge  released  by  the  shock  wave.  Solving  these  equations  gives 


qc(t)  = q(0  , t > 0 


(35) 


and 


qL(0 


C, 


c + c, 


q(t)  , t > 0 . 


(36) 


The  voltage  on  the  PZT  unit  and  load,  and  the  current  in  the  circuit  are, 
respectively, 


V(t)  = 


q(t) 
C + C, 


t 3*  0 


(37) 


and 


dqL(i)  CL  dq(t) 
i(t)  = — - = ~ — — , t > 0. 


dt 


C + CL  dt 


(38) 


Equations  (35)  through  (38)  can  be  solved  for  the  rectangular  or  hollow-cylinder 
unit  by  substituting  q(t)  from  Equations  (1)  or  (5).  respectively. 


B.  RECTANGULAR  PZT  UNIT 


Using  q(t)  from  Equation  (1)  in  Equations  (35)  through  (37)  gives 


I 


P A — - — (-)  , 0 < t 

0 C + CL  \t  / 


c 

P„A 

C + cL 


t > X . 


qL(t)  = 


( P‘'A  (x1)  ’ ° 

) P„A 

l ° ^'  + CL 


< t < r 


t > T 


V(t)  = 


I1)' 


0 < t < x 


t > X 


Figure  8(a)  is  a plot  of  the  normalized  charge  or  voltage  on  the  PZT  unit.  The 


charging  region  (0  < t < x)  is  characterized  by  a constant  slope  determined  by 
C',  /C.  For  the  open-circuit  case  (CL/C  = 0),  no  charge  is  transferred  to  the  load: 
this  is  identical  to  the  resistive  case  for  R C/x  = °°  as  shown  in  Figure  2(a).  For  the 
short-circuit  case  (C'L/C  = °°).  all  the  charge  is  transferred  to  the  load  capacitance  at 
zero  voltage:  this  situation  is  again  identical  to  the  resistive  case  lor  R C/x  = 0. 


The  circuit  current,  obtained  from  Equation  (38).  is 


iMiL.  . 0 < t < x 

’ x (’  + (’, 


t > X . 
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Figure  8.  Electrical  relations  for  a rectangular  PZT  unit  with  a capacitive  load  for 
various  values  of  CL/C.  (a)  Normalized  charge  qn(t)/P0A  or  voltage  V(t)/(PQA/C)  on 
the  PZT  unit,  (b)  Normalized  circuit  current  i(t)/(PQ A/r). 
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Figure  8(b)  is  a plot  of  the  normalized  current  for  the  rectangular  PZT  unit.  The 
current  is  a flat-topped  pulse  for  0 < t < r and  zero  for  t > t.  For  the 
short-circuit  case,  the  rate  of  charge  release  by  the  shock  wave,  P()A/t,  equals  the 
amplitude  of  the  short-circuit  current  pulse.  Fquation  (42)  and  Figure  8(b)  show 
that  the  amplitude  of  the  current  pulse  decreases  as  the  ('L/C  ratio  decreases  until, 
for  the  open-circuit  case,  the  current  is  zero. 

In  recent  experiments,  rectangular  PZT  56/44  and  PZT  95/5  units  were  depoled 
in  the  normal  mode  with  a capacitive  load.1  It  was  shown  that  tor  some  stress 
values  for  each  material  the  measured  voltage  and  current  pulses  were  similar  to 
those  in  Figure  (8). 

The  energy  on  the  PZT  unit  is  given  by 


lc-(t)  = 


j “swer 

r (PpAi-V  c y 
2f  \ c + c.  ) 


, 0 < t < T . 


t > r . 


and  the  energy  on  the  capacitive  load  is 


p;L(t)  = 


<1*0  A)-  C(|  /t 

2C  (C  + C,  ):  \t 


( P„  A)~  C C , 

2C  (C  + CL  )' 


, 0 ^ t < r 


t > r 


The  total  circuit  energy  I ( ( t ) = l ( (t)+l  ((t)  is 


( O’oA)2  / C' 


ft(.)=J  2c  vc  + c>. 


0 < t < r 


t > T 


Figures  9(a),  9(b),  and  9(c)  are,  respectively,  the  normalized  energy  on  the  PZT 
unit,  normalized  energy  on  the  load  capacitance,  and  the  total  normalized  energy. 
For  the  open-circuit  case,  all  the  released  charge  remains  on  the  PZT  unit;  the 
circuit  energy  has  a maximum  value  of  (P(1  A)2/2C.  For  the  short-circuit  ease  all  tile- 
energy  equations  are  zero  even  though  a current  P((A/r  Hows  in  the  circuit.  This  is 
analogous  to  the  resistive  load  tor  RC/r  = 0.  As  shown  in  Figure  9(b).  the  energy 
transferred  to  the  load  capacitance  is  a maximum  for  CL/C  = 1;  this  maximum 
value  is  (P()A)2/8C.  In  summary,  the  maximum  charge  is  transferred  to  the  load  lor 
C./C  = »,  but  the  maximum  energy  transfer  occurs  tor  CL/C  = 1. 


C.  HOLLOW-CYLINDER  PZT  UNIT 

In  this  section,  expressions  are  derived  tor  the  charge,  voltage,  current,  and 
energy  in  a circuit  consisting  ol  a hollow-cylinder  PZT  unit  and  a capacitive  load. 
The  charge  q ,(t>  on  the  PZT  unit  and  qL(t)  on  the  load  capacitance  are  obtained 
by  substituting  q(t)  from  liquation  (5)  into  liquations  (35)  and  (36).  respectively 
The  result  is 


p A — — — V-^— IT  . 0 < . < r 

■ fr-r 
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Figure  9.  Energy  relations  for  a rectangular  PZT  unit  with  a capacitive  load  for 
various  values  of  CL/C.  (a)  Normalized  energy  Ec  (t)/((P0 A)2  /2C)  on  the  PZT  unit, 
(b)  Normalized  energy  EL  (t)/((PQ  A)2/2C)  on  the  capacitive  load.  (c)  Total 
normalized  energy  E (t)/((P  A)2/2C). 
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Using  Equation  (37),  the  voltage  on  the  PZT  unit  and  load  capacitance  is 


Pn  A 


2ro  + t 


V(t)  = 


c+cL  in o +I 


T/  ^ t 
T 


0 < t < T, 


(48) 


PqA 

c + c, 


t > r . 


Figure  10(a)  is  a plot  of  the  normalized  charge  and  voltage  for  rQ/P  = 1.  In 
the  charging  region  (0  < t < r),  the  only  difference  between  these  curves  and  those 
for  the  rectangular  unit  in  Figure  8(a)  is  the  geometrical  factor 
[(2r0/fi)  + ( t/r) | /[(2r0/C)  + 1 ] . As  r() /t  the  factor  becomes  1 and  the  curves 

become  identical.  For  the  region  t > r,  the  curves  for  the  two  geometries  are  the 
same. 

The  current  equation  is 
PqA 


0 , t > t . 


Figure  10(b)  is  a plot  of  the  normalized  current  for  r()/C  = 1.  Again,  these  curves 
become  identical  to  those  for  the  rectangular  unit  in  Figure  8(b)  as  rQ/C  -*•«>.  At 
t/r  = 0.5,  the  normalized  current  for  the  hollow-cylinder  and  rectangular  units  are 
equal  for  any  value  of  r()/C  since  for  this  case  the  geometrical  factor  is  1. 

The  energy  on  the  PZT  unit  and  load  capacitance  are  given  by 


(b) 

Figure  10.  Electrical  relations  for  a hollow-cylinder  PZT  unit  with  a capacitive  load 
for  various  values  of  CL  1C.  Curves  are  given  for  rQ/l!  = 1.  (a)  Normalized  charge 
qc(t)/P0A  or  voltage  V(t)/(PQ A/C)  on  the  PZT  unit,  (b)  Normalized  circuit  current 
i(t)/(P0A/r). 


27 


and 


Figure  1 1 contains  the  normalized  energy  equations  for  the  hollow-cylinder  unit  tor 
r jSL  = 1.  The  only  difference  between  the  energy  equations  for  the  hollow-cylinder 
unit  (Equations  (50)  through  (52))  and  those  for  the  rectangular  unit 
(Equations  (43)  through  (45))  is  the  square  of  the  geometrical  factor. 


V.  INDUCTIVE  LOAD 


A.  GENERAL  CIRCUIT  EQUATIONS 

In  this  section,  the  load  in  Figures  1(a)  and  1(b)  is  an  inductance  L.  In  this 
case,  the  charge  will  oscillate  between  the  source  capacitance  and  load  inductance 
for  all  time  t.  Experiments  with  inductive  loads  have  shown  that  dielectric  relaxation 
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Figure  11.  Energy  relations  for  a hollow-cylinder  PZT  unit  with  a capacitive  load 
for  various  values  of  CL/C.  Curves  are  given  for  rQA!  = 1.  (a)  Normalized  energy 
Ec  (t)/((PQ  A)2 /2C)  on  the  PZT  unit,  (b)  Normalized  energy  EL  (t)/( (PQ  A)2 /2C)  on 
the  capacitive  load,  (c)  Total  normalized  energy  ET(t)/((P  A)2/2C). 


in  PZT  95/5  causes  severe  damping  of  the  circuit  current.5  Nevertheless  an  analysis 
is  presented  tor  the  idealized  case  for  completeness.  Experiments  with  other  PZT 
materials  may  be  in  closer  agreement  with  the  simple  theory. 

The  voltage  balance  equation  is 


q,(t)  _ di(t) 
C " L"dt 


t > 0. 


(53) 


Substituting  di(t)/dt  from  equation  (10)  gives  the  differential  equation 


d2qc(t)  qc(t)  _ d~ q( t ) 

dt2  LC  " dt2 


(54) 


This  equation  is  solved  for  q,(t)  by  substituting  q(t)  from  Equation  (1)  or  (5). 
Also  from  Equation  (10),  the  current  is 


i(t) 


dq(t)  dqc(t) 

dt  dt 


(55) 


B.  RECTANGULAR  PZT  UNIT 

Using  Equation  (1)  for  q(t)  gives  d2q(t)/dt2  = 0 for  all  t.  Equation  (54)  then 
becomes 


d2q,(t)  q,  ( t ) 
dt2  + LC 


t > 0 . 


(56) 


This  equation  is  first  solved  for  the  range  0 < t *5  r using  the  initial  conditions 
q.(0)  = i(0)  = 0.  From  Equation  (55).  i(t)  is  then  obtained.  To  find  q (t)  and  i(t) 
for  t > r.  a general  solution  is  obtained  for  q (t)  from  Equation  (56).  then  the 
current  is  obtained  from  the  equation  i(t)  = -dq  (t)'dt  (Since  dq(t)/dt  = 0 for 
t>r).  The  continuity  of  charge  and  current  at  t = r is  used  to  obtain  the  two 
constants  in  the  general  solution.  This  procedure  gives 


P0A  . 

sin  cot 


q,.(t) 


COT 

PqA 

COT 


0 < t < r . 

(sin  cot  - sin(cot  - co7)|  . t > t 


(57) 


PnA 


i(t) 


( 1 - cos  cot) 


0 < t < r . 


' T 


(58) 


[cos  (cot  - cot)  - cos  cot]  , t > T . 


In  these  equations  co  is  the  oscillatory  frequency  given  by  co2  = (LC)  ' . Using 
trigonometric  identities  to  simplify  the  q (t)  and  i(t)  expressions  gives 


PnA  / OOt  \ ZCOtX 

AsHtMt)  ■ 0 


^ t < T . 


q..(t) 


/ ^ sin  (-)  cohort-  t > 7, 

((f)  [:J  [ ^ 


(59) 


arid 


i(t) 


-Po  A 


-Po  A 


hr  . 0 < t < t . 

/U>T\  / COT\ 

in  ( — ) sinl  cot  - — 1 . t > t . 


(60) 


The  voltage  on  the  PZT  unit  is  obtained  from  the  equation  V(t)  = qc(t)/C.  Plots  of 
the  normalized  charge  and  voltage  are  given  in  Figure  12(a):  the  normalized  current 
is  given  in  Figure  12(b).  T = 27r/co  is  the  period  of  oscillation.  For  the  open-circuit 
case,  the  oscillatory  frequency  is  zero  (cot  = 0)  and  the  period  of  oscillation  is 
infinite  (T/t  = °°).  For  this  case  q(t)  = q(t)  and  i(t)  = 0 for  all  t.  This  can  be 
shown  for  cot  -*  0 by  using  the  relationships  sin(cot/2)/coT  -*■  t/2r.  cos(cot/2)  -*■  1 . and 
costcot  - cot/2  ) -*■  I in  Equation  (59)  for  q (t),  and  sin(cot/2)/r  -*  0 in  Equation  (60) 
for  i(t).  As  the  oscillatory  frequency  increases,  the  vibrational  amplitude  for  the 
normalized  charge  and  voltage  decreases.  For  the  short-circuit  condition  (cot  = °°). 
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Figure  12.  Electrical  relations 
various  values  of  cor  or 
V(t)/(PQ A/C)  on  the  PZT  unit. 


for  a rectangular  PZT  unit  with  an  inductive  load  for 
T It.  (a)  Normalized  charge  qc (t)/PQ A or  voltage 
(b) 


Normalized  circuit  current  i(t)/(PQ A/r). 
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the  frequency  of  oscillation  approaches  infinity,  but  the  amplitude  of  oscillation 
approaches  zero.  For  this  case,  it  can  be  shown  that  qc(t)  = 0 for  all  t by  using 
the  relationships  sin(cot/2)/co7  -*■  0 and  sin(cor/2)/<or  -*■  0 for  uit in 
Equation  (59).  For  0 < t < r the  current  oscillates  about  P(JA It  with  amplitude 
P A/r:  for  t > t the  current  oscillates  about  a zero  value  with  a variable  amplitude 
(2P„  A/t)  sin(cor/2).  When  c or  is  an  odd  multiple  of  7r,  the  amplitude  is  2P0A/r. 
When  cor  is  an  even  multiple  of  n.  the  amplitude  is  zero.  The  short-circuit  current 
(cor  = °°)  can  be  obtained  from  Equation  (55)  using  qc(t)  = 0 for  all  t,  or  from 
Equation  (60)  using  sin’(<ot/2)  -*  1/2  and  sin(cor/2)  -*■  0 as  cor->°°.  As  expected  the 
result  is 


i(t.COT  = oo)  = 


i P„A 
1 r 

I o 


0 < t . < r , 

t > r . 


(61) 


Expressions  for  the  energy  on  the  PZT  unit  and  load  inductance  can  be 
obtained,  respectively,  from  the  equations  tv(t)  = q^(t)/2C  and  EL(t)  = L i 2 ( t )/ 2 . 


The  results  are 
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and 
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The  expression  for  the  total  circuit  energy  ET(t)  = F(.  (t ) + E,  (t)  is 


Figures  13(a),  13(b),  and  13(c)  are  plots  of  the  normalized  energy  expressions 
for  various  values  of  the  parameter  c or.  For  the  open-circuit  condition  (cor  = 0).  it  can  be 
shown  that  FT't)  = [ ( P(> A)2 /2C] (t/r)2  for  0 < t < r and  ET(t)  = (P()A)2/2C  for 
t > r.  and  that  E,  (t)  = 0 for  all  t.  Also,  for  the  short-circuit  condition  (cor  = °°). 
Ec(t)  = EL(t)  = HT(t)  for  all  t.  Figures  13(a).  13(b).  and  13(c)  indicate  that  for 
0 < t < r the  PZT.  inductive,  and  total  energies  are  oscillatory  in  nature.  For 

t > r,  the  total  circuit  energy  is  constant,  but  it  oscillates  between  the  PZT  unit 
and  inductance  with  frequency  co.  This  situation  can  be  compared  with  the 
capacitive  load  in  which  the  energies  on  the  PZT  unit  and  load  capacitance 
remained  Fixed  for  t > r.  The  total  energy  curves  for  cor  = 0 and  cor  = 0. 1 are 
within  about  0.19?  of  each  other  and,  therefore,  could  not  be  shown  as  separated 
on  the  plot. 

For  t = r and  cor  = 1 (Figure  13(a)).  the  normalized  PZT  energy  is  about 
70'/  of  the  maximum  theoretical  value.  At  t = 2.1r.  the  PZT  energy  is  zero 

(corresponding  to  the  case  of  zero  charge  and  voltage  on  the  PZT  unit  and 

maximum  circuit  current  (see  Figures  12(a)  and  12(b)).  but  at  about  t = 2.h5r  the 
PZT  energy  has  increased  to  about  90'/.  since,  at  this  time  the  energy  that  was  on 
the  inductor  at  t = r (about  20' V' ) has  been  transferred  back  to  the  PZT  unit.  For 
this  case,  considerable  charge  (about  0.95  P(| A)  is  being  transferred  back  and  forth 
between  the  PZT  unit  and  inductance  (Figure  12(a)).  The  inductive  energy  is  about 
90'?  of  the  theoretical  maximum  value  at  t = 2. 1 r . (This  corresponds  to  the  case  of 
zero  capacitive  energy.)  This  inductive  energy  transfer  is  greater  than  the  25'? 

maximum  energy  transfer  for  the  capacitive  load  (Figure  9(b)). 

Note  that  in  Figure  13(b)  the  wt  = I curve  peaked  at  t = 2.lr.  As  tor 
decreases  from  (his  value,  the  amplitude  of  the  curve  increases,  and  the  peak  occurs 
for  larger  times.  For  cor  = 0.1.  the  amplitude  of  the  peak  is  99';  of  the  theoretical 
maximum  value  and  occurs  at  t = 164r.  As  cor  increases  from  cor  = 1.  the 
amplitude  of  the  peak  decreases  and  occurs  for  smaller  times.  The  result  for 
cor  = 10  is  shown  in  the  figure. 
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Figure  13.  Energy  relations  for  a rectangular  PZT  unit  with  an  inductive  load  for 
various  values  of  tor  or  T It.  (a)  Normalized  energy  Ec  (t)/((PQ  A)Z/2C)  on  the  PZT 
unit,  (b)  Normalized  energy  EL  (t)/((PQ  A)2/2C)  on  the  inductive  load,  (c)  Total 
normalized  energy  ET  (t)/( (PQ  A)2 /2C). 


C.  HOLLOW-CYLINDER  PZT  UNIT 


Substituting  d q(t),dt"  Irom  Equation  (5)  into  Equation  (54)  gives  the 
differential  equations  for  q (t). 


d^L(t)  q^O 
dt2  LC 


0 < t < r 


d2qc(t)  qL.(t) 
dt  + LC  “ ° 


t > T 


The  same  procedure  is  used  lor  solving  these  equations  ami  Equation  (55)  for  q (t) 

and  i(t)  that  was  used  for  solving  the  equations  for  the  rectangular  IV.T  unit.  This 
procedure  gives 
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Using  trigonometric  identities  gives 
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Plots  of  the  normalized  charge  and  voltage  are  given  in  Figure  14(a)  for  r(t/C  = 1; 
the  normalized  charge  is  given  in  Figure  14(b). 
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Equations  (69)  and  (70)  arc  both  continuous  at  t = t.  Also,  these  equations 
reduce  to  Equations  (59)  and  (60),  respectively,  lor  the  rectangular  PZT  unit  as 
r()/e->-°°.  Since  (2r0/C)/[(2r0/C)  + 1 1 ->•  1,  and  l/(rn/fi)  -*■  0,  as  r 0/e->°°  only  the  first 
terms  in  Equations  (69)  and  (70)  are  nonzero,  and  these  are  just  the  terms  in 
Equations  (59)  and  (60)  for  the  rectangular  PZT  unit. 

For  the  open-circuit  case  (cor  = 0),  the  limiting  values  as  cor  -*■  0 for  the 
trigonometric  relationships  given  in  the  previous  section  for  the  rectangular  unit  can 
be  used  to  show  that  Equations  (69)  and  (70)  become,  respectively. 


0 < t < r , 

(71) 

t > r , 
and 


i(t,cor  = 0)  = 0 , t 3s  0 . 


(72) 


Equation  (71)  is  just  the  charge  release  q(t)  from  a hollow-cylinder  PZT  unit 
(Equation  (5)). 


For  the  short-circuit  case  (cor  = °°),  the  relationships  sin(cot/2)  -*•  0. 
sin2 (cot/2) -»  1/2  as  cor-*'00  can  be  used  to  show  that  Equations  (69)  and  (70) 
reduce  to 


q.(t.cor  = °°)  = 0 . t > 0 


(73) 


and 


( [>„A  V g r 

) - (iLo 


i(t.cor  = °°) 


+ 1 


C < t < r 


(74) 


0 


t > r . 


As  expected,  the  short-circuit  current  in  Equation  (74)  agrees  with  the  short-circuit 
current  for  the  resistive  load  (Equation  (27))  and  capacitive  load  (Equation  (49) 
with  CL/C  = °°).  The  voltage  on  the  hollow-cylinder  unit  can  be  obtained  by 
substituting  q c ( t ) from  Equation  (69)  into  the  equation  V(t)  = q,(t)/C. 

Expressions  for  the  energy  on  the  PZT  unit  and  inductance  are  obtained, 
respectively,  from  Ec(t)  = q^(t)/2C  and  EL(t)  = L i2 (t)/2 . Making  substitutions  for 
qc(t)  from  Equation  (69)  and  i(t)  from  Equation  (70)  gives 
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The  total  circuit  energy  is  obtained  by  adding  Equations  (75)  and  (76),  to  give 
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These  energy  expressions  are  continuous  at  t = r.  Also,  as  r()/8-*°°  Equations  (75), 
(76),  and  (77)  reduce  to  Equations  (62),  (63),  and  (64)  for  the  rectangular  PZT 
unit,  respectively.  Figures  15(a),  15(b),  and  15(c)  are  plots  of  the  normalized  energy 
expressions  for  r0/C  = 1. 


VE  GEOMETRICAL  SIZE  ESTIMATES  OF  PZT  UNITS  IN 
FERROELECTRIC  POWER  SUPPLIES  FOR  SOME  EXAMPLE 
HIGH-VOLTAGE  ENERGY  CONVERTERS 


In  this  section,  the  geometrical  size  of  PZT  units  for  use  in  pulsed  power 
supplies  will  be  estimated  for  some  example  energy  converters.  Energy  converters 
with  resistive  and  capacitive  input  impedances  were  chosen.  The  electrical  equations 
that  were  derived  in  Sections  111  and  IV  will  be  used  in  the  analysis.  Expressions 
for  the  area  A and  height  h of  the  PZT  unit  will  be  obtained  in  terms  of  the 
converter  input  resistance  R or  capacitance  C'L  and  the  input  voltage  V.  The  shock 
and  electrical  parameters  for  the  PZT  unit  will  be  taken  as  those  for  PZT  56/44 
material.  A 12.6-mm  cube  of  PZT  56/44  had  a capacitance  C of  240  pF.1  For  this 
cube,  it  was  shown  experimentally  that  about  40  pC  of  the  51  pC  of  available 
charge  was  released  at  a maximum  electric  field  F(|  of  about  2 kV/mm  in  the  7 to 
8 GPa  stress  range.1  In  this  analysis,  P()  will  therefore  be  taken  as  0.25  fiC/mnr 
(40  pC/(12.6  mm)2).  (The  PZT  56/44  remanent  polarization  vector  PQ  is 
0.33  juC/mm2.1  We  have  simplified  the  present  analysis  by  taking  P0  as  representing 
the  actual  charge  released.)  The  average  measured  shock  velocity  U was  4.65  km/s.1 
The  geometrical  shape  of  the  PZT  unit  will  be  discussed  after  the  size  has  been 
estimated.  It  should  be  noted  that  in  most  experiments,  stresses  from  relief  wave 
interactions  will  have  destroyed  the  unit  for  times  significantly  greater  than  the 
shock  transit  time  r. 


A.  CONVERTERS  WITH  RESISTIVE  INPUT  IMPEDANCES 

Table  1 lists  two  converters  and  their  characteristics.  These  converters  require 
approximately  fiat-topped  input  voltage  pulses.  Figures  2(a)  and  4(a)  for  the 
rectangular  and  hollow-cylinder  units,  respectively,  indicate  that  the  value  of  RC/r 
must  be  near  that  for  short-circuit  conditions  (0.1  or  less)  for  this  type  of  pulse. 
The  load  voltage  reaches  a constant  value  soon  after  the  shock  wave  enters  the  PZT 
unit  and  does  not  begin  to  decrease  until  after  the  shock  wave  passes  out  of  the 
unit.  Expressions  for  the  voltage  near  short-circuit  conditions  (RC/t -*•())  can  be 
obtained  using  q.(t)  from  Equation  (12)  for  the  rectangular  unit  and  using  qc(t) 
from  Equation  (23)  for  the  hollow-cylinder  unit  in  the  equation  V(t)  = q.(t)/C.  The 
results  are 
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Figure  15.  Energy  expressions  for  a hollow-cylinder  PZT  unit  with  an  inductive  load 
for  various  values  of  cor  or  T It.  Curves  are  given  for  rQ/t  =1.  (a)  Normalized 
energy  Ec  (t)/((PQ  A)2/2C)  on  the  PZT  unit,  (b)  Normalized  energy  EL  (t)/((PQ  A)2/2C) 
on  the  inductive  load,  (c)  Total  normalized  energy  ET  (t)/((PQ  A)2 /2C). 
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Table  1.  Energy  converters  with  resistive  input 

impedances. 1 3 

Pulse  Duration  of 

Input 

Input 

Flat-Topped 

Resistance  R 

Voltage  V 

Input  Pulse  r 

Converter 

(kF2l 

(kV) 

(/as) 

Magnetron 

0.4  - 1 

50 

0.1  - 10  (use  2.7  /as)a 

Klystron 

0.4  - 2 

100  - 1000 

(use  2.7  /as)a 

■'T 'ho  2.7-  n'  value  lor  the  pulsi  duration  is  Ills-  average  transit  time  r calculated  using  a 4.65  km/s  shock  transit  time  for 
IV  I 56/44  material  and  a PZT  unit  thickness  S of  1 2.6  mm. 


lor  the  hollow-cylinder  unit.  For  r 0/C  > 1.  Hquation  (79)  reduces  to  Equation  (78). 
Equation  (78)  will  be  used  for  estimating  the  size  of  the  PZT  unit.  Figure  2(a)  for 
RC/r  = 0.1  indicates  that  the  normalized  charge  and  voltage  are  about  1 09?  of  their 
maximum  values.  This  is  also  indicated  in  Equation  (78)  by  substituting  0.1  for 
RC/r. 

The  height  h.  thickness  C.  and  width  w for  the  PZT  unit  can  be  calculated 
using  Equation  (78)  for  V and  the  electrical  and  shock  characteristics  of  the  PZT 
56/44  material. 


r- 


i iii 


The  height  h is  given  by 


h = V/F0 


where  FQ  = 2 kV/mm  and  V is  the  applied  voltage.  The  thickness  C is  given  by 


8 = Ur 


where  U = 4.65  km/s  and  r is  the  required  input  pulse  duration.  The  width  w 
obtained  from  Equations  (78)  and  (81)  is 


RP0U 


where  R is  the  converter  input  impedance  and  P() 
the  capacitance  of  a 12.6-mm  PZT  56/44  cube,  then 


= 0.25  nC/mm2.  If  240  pF  is 


C = 0.019- 
h 


is  the  capacitance  in  nF  of  a PZT  56/44  unit  witli  area  A in  mnr  and  height  h in 


The  dimensionless  parameter  RC/r  in  Equation  (78)  can  be  obtained  using  V from 
Equation  (80)  and  C from  Equation  (83).  The  result  is 


R(  F„ 

— = 0.0  9—  = 0.15 
r P„ 


where  the  values  for  F()  and  Pn  have  been  used.  For  this  value  of  RC/r  about  15% 
of  the  available  charge  will  remain  on  the  PZT  unit  and  85%  flow  through  the 
resistive  load  (Figure  2(a)).  At  t = r,  the  energy  EL(r)  dissipated  in  the  load 
resistance  R is  given  by 


44 


I'll  is  can  be  shown  using  i = V/R  in  the  equation  L£l  ( t ) - I i'(t)dt. 

Jo 

In  the  analysis  to  follow,  the  thickness  8 of  the  PZT  unit  will  be  taken  to  be 
12.6  mm  giving  a pulse  duration  r of  2.7  ps.  This  thickness  was  chosen  because  it 
corresponds  to  a measured  charge  release  of  about  77%. 1 If  a different  thickness 
were  used,  then  the  charge  release  may  be  different,  requiring  a PQ  different  from 
Q.25  pC/ mm'  in  the  calculations. 

Table  2 gives  the  PZT  geometrical  size  estimates  for  the  energy  converters.  The 
input  resistances  and  voltages  were  chosen  to  cover  the  range  of  input  characteristics 
in  Table  1.  The  size  estimates  for  the  magnetron  indicate  that  a single  rectangular 
unit  or  a series  of  rectangular  units  connected  electrically  in  parallel  would  provide 
a reasonable  power  source.  In  Reference  1,  a shock-depoled  PZT  56/44  unit 
25.9-mm  high,  25. 6-nnn  wide,  and  12.4-mm  long  (unit  volume  8 cm3)  was  used  to 
pulse  charge  a 1.2  nF  load  capacitor  to  55  kV.  About  78%  of  the  available  charge 
was  released.  This  unit  is  similar  in  size  to  the  estimated  PZT  units  for  the 
magnetron  in  Table  2.  The  size  estimates  for  the  klystron  indicate  that 
hollow-cylinder  units  may  be  more  appropriate  than  rectangular  units  due  to  the 
large  amount  of  PZT  material  required.  Also,  the  large  height  of  PZT  material 
(500  mm)  for  the  1000-kV  input  voltage  may  be  prohibitive  due  to  the  extreme 
care  that  would  be  required  during  fabrication  of  the  unit  from  the  PZT  elements. 

Note  that  in  the  calculations  for  Table  2 the  height  h and  width  w are 
determined  from  the  input  resistance  and  voltage,  but  the  length  8 depends  only  on 
the  pulse  duration  t.  A change  in  t affects  the  volume  of  PZT  material  required  to 
drive  the  converters,  but  not  the  shape  of  the  voltage  pulse  since  RC/t  is  constant. 
Also,  as  expected,  changing  r also  changes  the  energy  dissipated  in  the  resistance 
(liquation  (85)). 


B.  CONVERTERS  WITH  CAPACITIVE  INPUT  IMPEDANCES 


Table  3 lists  five  converters  and  their  characteristics.  Expressions  for  the  area  A 
and  height  h of  the  PZT  unit  will  be  calculated  in  terms  of  the  converter  input 
capacitance  C,  and  input  voltage  V.  The  voltage  on  the  PZT  unit  and  converter  for 
t = r will  be  used  for  estimating  the  size  of  the  unit.  Equations  (41)  and  (48)  are 


Table  2.  Estimated  geometrical  and  electrical  characteristics  of  PZT 
units  for  energy  converters  with  resistive  input  impedances. 


the  expressions  for  the  rectangular  and  hoi  low -cylinder  units,  respectively.  At  t - r, 
the  voltage  expressions  for  the  two  geometries  are  the  same.  The  expression  is 


Po  A 

t = T . 

capacitive  input  impedances.1 3 

(86) 

Table 

V - , 

c + cL 

3.  Energy  converters  with  i 

Input 

Input 

Operating 

Capacitance  CL 

Voltage  V 

Frequency 

Converter 

InF) 

(kV) 

(MHz) 

Pancake  ( LC) 

1 - 300 

500  - 5000 

0.5  - 50 

Oscillator 

Spiral 

1 000  - 1 0000 

10-50 

— 

Oscillator 

Frozen-Wave 

0.1-10 

20-100 

10-  1000 

Generator 

Ross  Circuit 

0.01  - 1 

20-100 

10-  16000 

Transmission 

0.01  - 1 

20  - 100 

50  - 1 000 

Line  Oscillator 

The  height  of  the  PZT  unit  is  given  by  Equation  (80).  The  capacitance  C of  the 
unit  is  given  by  Equation  (83).  Substituting  Equation  (80)  into  Equation  (83)  gives 


CV  = 0.038A  (87) 


% 


47 


with  C in  nF,  V in  kV,  CV  in  j uC,  and  A in  mm2.  Substituting  Equation  (87)  into 
Equation  (86)  gives  an  expression  for  the  area  A.  The  result  is 


A = 4.72  CLV 


(88) 


with  CL  in  nF  and  CLV  in  pC.  Using  the  value  for  CL  and  V from  Table  3.  the 
size  of  the  PZT  unit  can  be  estimated  from  Equations  (80)  and  (88).  The  ratio 
C L/C  can  be  obtained  by  dividing  Equation  (87)  by  Equation  (88).  The  result  is 


Cl 

C 


= 5.56  . 


(89) 


Table  4 gives  the  PZT  geometrical  size  estimates  for  the  capacitive  energy 
converters.  The  thickness  C was  taken  as  12.6  mm  with  a corresponding  pulse 

duration  of  2.7  ps  (as  for  the  converters  with  resistive  input  impedances).  The  input 

capacitances  and  voltages  cover  the  range  of  values  listed  in  Table  3.  Since  r()/£  > 1 
both  the  rectangular  and  hollow-cylinder  units  give  the  same  response  for  all  the 
cases  shown.  Table  4 indicates  that  the  PZT  units  for  the  maximum  input 
characteristics  for  the  pancake  and  spiral  oscillators  are  unreasonably  large  for  power 
supply  fabrication.  It  may  be  possible  to  fabricate  the  required  PZT  power  supplies 

for  the  other  listed  input  characteristics,  but  in  some  cases  it  may  be  impractical. 

For  the  spiral  oscillator,  with  input  characteristics  of  1000  nF  and  10  kV,  it  may 
be  possible  to  fabricate  a power  supply  by  stacking  a series  of  hollow-cylinder  units 
(each  5 mm  in  height)  that  are  electrically  in  parallel,  each  separated  with  a thin 
layer  of  insulating  material.  The  0.75-mm  width  ot  the  PZT  unit  lor  the 
frozen-wave  generator  and  Ross  circuit  with  input  characteristics  of  0.1  nF  and 
20  kV  is  probably  impractical  and  unnecessary.  If  the  input  voltage  were  increased 
to  100  kV  the  width  would  increase  to  a more  reasonable  value  of  3.7  mm.  For 
both  the  resistive  and  capacitive  input  impedance  converters,  it  may  also  be  possible 
to  achieve  a multishot  power  supply  capability  via  the  sequential  depoling  of 
multiple  units. 


Table  4.  Estimated  geometrical  and  electrical  characteristics  of  PZT  units  for 
energy  converters  with  capacitive  input  impedances. 


Vi.  SUMMARY 


Electrical  equations  have  been  derived  for  the  normal-mode  shock-depoling  of 
PZT  ferroelectric  ceramics  with  resistive,  capacitive,  and  inductive  loads.  To  simplify 
the  analysis,  it  was  assumed  that  the  PZ.T  material  responded  in  an  ideal  manner; 

i.e.,  charge  was  released  completely  and  instantaneously  at  the  shock  front. 
Rectangular  and  hollow-cylinder  shapes  were  considered  for  the  PZT  unit.  Electrical 
expressions  were  obtained  for  the  charge,  voltage,  current,  and  energy  in  the  circuit 
as  the  shock  front  passed  through  the  PZT  unit  and  after  the  shock  front  passed 
out  of  the  unit.  Plots  of  these  equations  were  presented  as  the  magnitude  of  the 
load  varied  from  short-circuit  conditions  to  open-circuit  conditions.  These  results  and 
the  electrical  properties  of  shock-depoled  PZT  56/44  ferroelectric  ceramics  were  used 
to  estimate  the  geometrical  size  and  electrical  characteristics  of  PZT  units  for 
selected  energy  converters  with  resistive  and  capacitive  input  impedances. 
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